Weak-Field Thermal Hall Conductivity in the Mixed State of rf-Wave Superconductors 
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Thermal transport in the mixed state of a d-wave superconductor is considered within the weak- 
field regime. We express the thermal conductivity, k xx , and the thermal Hall conductivity, K xy , in 
terms of the cross section for quasiparticle scattering from a single vortex. Solving for the cross 
section (neglecting the Berry phase contribution and the anisotropy of the gap nodes), we obtain 
k X x(H,T) and K xy (H,T) in surprisingly good agreement with the qualitative features of the experi- 
mental results for YBa2Cu3 06.99- In particular, we show that the simple, yet previously unexpected, 
weak-field behavior, K xy (H,T) ~ TyH, is that of thermally-excited nodal quasiparticles, scattering 
primarily from impurities, with a small skew component provided by vortex scattering. 

PACS numbers: 74.25.Fy, 74.60.Ec, 74.72.-h 



Thermal Hall conductivity provides the most direct 
measure of low temperature quasiparticle transport in 
a d-wave superconductor. Since quasiparticles are part 
electron and part hole, their energy is well defined but 
their charge is not. Thus, it is thermal current that fol- 
lows quasiparticle current. The longitudinal thermal con- 
ductivity, k xx , has both an electronic and a phononic con- 
tribution. However, the thermal Hall conductivity, K xy , 
induced by a perpendicular magnetic field (the Righi- 
Leduc effect), is purely electronic in origin and the direct 
consequence of a transverse quasiparticle current. 

Over the past few years, much progress has been made 
in measuring the thermal Hall conductivity of the cuprate 
superconductors in the mixed (vortex) state, H c \ < H < 
H C 2 0, Hill Most recently, Ong and co-workers H 
measured n xy in high-purity single crystals of slightly 
overdoped (T c = 89 K) YBa 2 Cu 3 6 .99 (YBCO). Their 
data indicates that, for magnetic fields up to 14 Tesla 
and temperatures between 15 K and 28 K, K. xy /T 2 is 



only a function of the ratio v H/T. This is in agreement 
with the scaling theory proposed by Simon and Lee 
which predicts that, for nodal quasiparticles with a Dirac- 
like dispersion, n xy (H, T) ~ T 2 F xy (VH/~fT) where 7 = 
(ks /vf)y/c/He and F xy (x) is a general scaling function. 
Furthermore, the experiments show that for y/H <C "fT 
the measured scaling function has the surprisingly simple 



form, F xy (x) ~ x, such that n xy (H,T) — CoT^/H where 
Co is a constant. For larger magnetic fields, the measured 
curves peak and then decrease. It is unusual to see \/H 
(rather than H) in a Hall response and this interesting 
result was theoretically unexpected. Yet such a simple 
functional form must have a simple explanation. In this 
Letter, we seek to provide it. 

The key to this problem is that we are dealing with the 
high temperature regime of low T quasiparticle transport. 
By low T quasiparticle transport, we mean that quasi- 
particles are excited only in the vicinity of gap nodes 
and inelastic scattering can be neglected. (Quasiparti- 
cle dispersion is therefore given by the anisotropic Dirac 
spectrum, E = (v 2 k 2 + wf^l) > where Vf is the Fermi 



velocity, V2 is the slope of the gap, and fej and ki are de- 
fined locally about each node. With our choice of axes, 
gap nodes are located at ±pi?x and ±pfY in momentum 
space.) However, the experiments involve temperatures 
large compared to the impurity scattering rate and the 
vortex scattering rate. As a result, the quasiparticles 
responsible for transport are thermally generated rather 
than impurity-induced 0, Q> @ or magnetic field- induced 
0, 0. (Note that thermal transport in the opposite, 
low T, regime has been discussed frequently in the re- 
cent literature [HI El El HE!.) This hi g h T regime is 
relatively simple. To understand the thermal conductiv- 
ity, we need only understand how the thermally excited 
quasiparticles scatter from impurities and magnetic vor- 
tices. 

In the presence of a magnetic field {H& < H <C H C 2), 
vortices penetrate the sample (a 2D Cu02 layer). They 
are distributed randomly, pinned to local defects. The 
cuprates are extreme type II superconductors in which 
the coherence length, £, is much smaller than the penetra- 
tion depth, A. As a result, while the vortex cores may be 
well separated, the magnetic field profiles overlap signif- 
icantly such that there is little variation in the magnetic 
field across the sample. We therefore adopt the extreme 
type II limit of £ — > and A — > 00 and take the magnetic 
field to be constant, H = Hi. In this limit, there are 
only two remaining length scales. The first, 1/fe, is set 
by the temperature such that k = E/hvj = krjT/hvf. 
The second length scale, R, is half of the average dis- 
tance between vortices. With one flux quantum per vor- 
tex, HttR 2 = $0 = hc/2e, so we define R = J hc/eH. 
In terms of R, we can define the (2D) density of vortices 
to be n v — H/&o — 1/ttR 2 . The ratio of the two length 
scales yields kR = yT/vH, which is the inverse of the 
argument of the scaling functions. 

We consider nodal quasiparticles carrying a heat cur- 
rent in response to a thermal gradient in the x-direction. 
Defining a mean free path, £, we can express k xx in terms 
of the electronic specific heat, C v , via k xx = v£C v /2 
where v is the average quasiparticle velocity |l6( ■ Specific 
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heat in the mixed state has been calculated by Kopnin 
and Volovik |l(J . For magnetic fields small compared to 
the temperature, they find C v ~ T 2 [l + 0{{V~H ' llT) 2 )]. 
The second term, a measure of the magnetic field con- 
tribution to quasiparticle generation, can be neglected 
in the regime of interest. Doing so and defining a 
thermal Hall angle, tan#y = k xv /k xx , yields a simple 
form for the thermal conductivity: k xx /T = a^kl and 
Kxy/T — aokitanOn, where ao ~ 1.72k%v/tiV2 times 
the stacking density of CuC>2 planes. 

The mean free path has contributions from both impu- 
rity scattering and vortex scattering. For small impurity 
densities and dilute vortices, we expect these to be rel- 
atively independent. Thus, via Matthiessen's rule, we 
write l/l = 1/£q + l/£ v where £q and £ v are the contri- 
butions from impurities and vortices respectively. Since 
l/£ v vanishes for H = 0, £q can be found empirically 
from K e (T), the electronic part of k xx (H = 0). We de- 
fine A = k£ — n e (T)/aoT. As argued by Simon and 
Lee , we expect scaling even in the presence of disor- 
der as long as the impurity scattering does not yield an 
additional length scale (as for ^-correlated disorder). If 
so, then A must be T- independent. (Note that it is a dif- 
ferent constant from that obtained in the universal limit 
0.) Zero-field measurements of k xx in YBCO 0, 0] 
show that this is realized experimentally for T < 30 K, 
which is precisely the temperature range over which the 
n xy data obeys scaling. 

Now consider the vortex contribution. Superflow cir- 
culates about the vortices, falling off like 1/r near the 
center of each vortex and cancelling with superflow of op- 
posite direction (circulating about neighboring vortices) 
in the regions between vortices. We therefore attribute 
to each vortex the circulating superflow in the area that 
surrounds it. Since the vortex scattering should depend 
primarily on the average vortex density and be relatively 
insensitive to higher order details of the vortex configu- 
ration, we can approximate the exact superflow distribu- 
tion by cutting off the superflow about each vortex at a 
distance R from its center. In this manner, we define an 
effective single vortex for which we have approximately 
accounted for the influence of neighboring vortices. If 
we assume that scattering events from such vortices are 
otherwise uncorrelated, then we can express £ v in terms 
of the single vortex transport cross section, a\\ , and the 
density of vortices, n v . In this approximation, the trans- 
port scattering rate is 1/t v = n v vau and, since £ v = vt v , 
we find £ v = l/n v an = irR 2 /a\\. Since vortices are en- 
dowed with a circulation, the vortex scattering cross sec- 
tion can have a small skew component, a±. Thus, when a 
quasiparticle encounters a vortex, it has a slightly greater 
probability of scattering to one side than the other and 
thereby contributes to J y . This process repeats with each 
successive vortex until the quasiparticle has travelled a 
distance equal to its mean free path. Therefore, given 
a heat current, J x , we can express the transverse heat 



current as J y — J x n v a±£. The thermal Hall angle is 
then given by tan## = n xy /n xx = —J y /J x = —<j±£/itR 2 
where we used n xy — ~n yx . 

Combining our results to this point, we find that k xx 
and K xy can be expressed in terms of the single vortex 
scattering cross section via 
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where x = 1/kR = y/H/jT, Ep = Vfpp/2, and we have 
defined dimensionless functions of kR = 1/x such that 
fc<T|| = /||(/ci?) and ka± = —(k/pp)f±(kR). The extra 
factor of k/pF in the definition of /j_ reflects the fact 
that K xy is small by a factor of ksT '/ ' Ep 0. 

Next, we calculate /||(fci?) and f±(kR) by consider- 
ing the quantum mechanical scattering of a quasiparti- 
cle from a single vortex. Our calculation, the details 
of which will be presented elsewhere 18j, is similar in 
spirit to that conducted by Cleary [ijj for the case of 
an s-wave superconductor. We consider the Bogoliubov- 
de Gennes (BdG) equation for a ci-wave superconductor 
in the presence of a single vortex. We apply a singular 
gauge transformation that simplifies the BdG Hamilto- 
nian at the cost of imposing antiperiodic boundary con- 
ditions which require that our wave function change sign 
with each trip around the vortex. We further simplify by 
shifting the origin of momentum space to the location of 
one of the gap nodes and neglecting scattering from one 
node to another. (This is physically reasonable since the 
superflow from which quasiparticles scatter is smooth on 
the scale of 1/pp.) The resulting problem is one of an 
(anisotropic) Dirac fermion scattering from an effective 
non-central potential (due to the superflow) in the pres- 
ence of antiperiodic boundary conditions and small, yet 
important, curvature terms in the Hamiltonian. 

Quasiparticles interact with vortices via the superflow 
as well as the Berry phase factor of (-1) acquired upon 
circling a vortex. This phase is encoded in the antiperi- 
odic boundary conditions imposed on quasiparticles in 
our chosen gauge. We make the following approxima- 
tions. First, we neglect the Berry phase effect and in- 
stead adopt periodic boundary conditions for the quasi- 
particle, which is equivalent to considering the case of an 
hc/e (double) vortex. To justify this, we note that the 
Berry phase only affects quasiparticle trajectories that 
lie within the thermal deBroglie wavelength of the vor- 
tex core, such that two paths that pass on cither side 
of the core can interfere. Thus, any additional contribu- 
tion to the cross section can at best equal the deBroglie 
wavelength. Since this wavelength is much smaller than 
the spacing between vortices (for \[H <C 7T), the Berry 
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phase contribution may change the size of the cross sec- 
tion, but should not change its dependence on magnetic 
field. Second, we assume an isotropic Dirac dispersion in 
the linearized Hamiltonian, which allows us to more eas- 
ily work with angular momentum eigenstates and simpli- 
fies the calculation. This is clearly an approximation for 
the cuprates where Vf exceeds V2 by a factor of 10 to 20 
|20j |. Third, we approximate the effect of neighboring vor- 
tices by cutting off the superflow distribution about our 
single vortex at a distance R from its center. By construc- 
tion, the flux through this circle is exactly one (hc/2e) 
flux quantum. The resulting superfluid momentum (su- 
perflow) is P s (r) = (h/2)(l/r-r/R 2 )6(R-r)4>. The BdG 
Hamiltonian (for quasiparticles about the node at +pf'X-) 
is then given by the sum of a linearized (Dirac) part, 
Hd = v f[ T 3Px + T\Py + P sx ], and a quadratic (curvature) 
part, H c = (v f /2 PF )[ T3 (p 2 + P 2 ) + 2P S ■ p + Ti(2p xPy )], 
which is small sufficiently far from the vortex center 
(r > 1/pf)- Since we cutoff our model at the scale of 
the vortex core (£ ~ 10/ P f), curvature terms can be 
considered perturbatively. Finally, we select a reason- 
able core size and model the vortex core as a region with 
vanishing superflow; which is the best we can do in the 
absence of further experimental input. We now have a 
well-defined scattering problem, which is first solved con- 
sidering the linearized Hamiltonian, and then perturbed 
to first order in the curvature terms. (Note, if curvature 
terms are neg lected completely, there is no skew scatter- 
ing jEElEll-) Due to these approximations, we cannot 
expect our results to be quantitative. Yet, the qualitative 
agreement with experiment is surprisingly good, which 
indicates that we have retained the essential physics. 

Note that Hr> includes an effective non-central poten- 
tial, V = —VfP s (r) sin cj>. This mixes angular momentum 
eigenstates and requires that we solve for all of them 
simultaneously. We do so numerically, including the con- 
tributions of up to 46 angular momenta. Given the eigen- 
states both inside (r < R) and outside (r > R) the vor- 
tex, we apply boundary conditions at the origin, match 
solutions at the vortex edge (r = R), and construct a 
wave function composed of an incident plane wave and 
an outgoing radial wave. The angular prefactor of the ra- 
dial piece yields the differential cross section. Summing 
over contributions from each of the four nodes and inte- 
grating over scattering angles, we obtain the total cross 
section, a = J dip (da / 'dip) , the transport cross section, 
eni = J dip (1 — cosy) (da /dip), and the skew cross sec- 
tion, a± = J dip sin ip (da / 'dip) . 

Results for a range of intervortex distances, kR = 
jT/VH, from 0.5 to 15, are plotted in Fig. □ Note 
that while a and a^ are plotted in units of 1/fc, a± is 
plotted in units of —1/pp- This reflects the fact that 
the skew cross section, induced by the curvature terms 
of the Hamiltonian, is small by a factor of k/pF- The 
minus sign indicates that the quasiparticles get deflected 
to the right, just as an electron would in response to the 
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FIG. 1: Calculated total, transport, and skew cross sections 
as a function of kR from 0.5 to 15. For a and cr±, solid lines 
are fits to straight lines with small negative intercepts. For 
<7|l , the solid line is an interpolation of the numerical data in 
which the constant large-fei? behavior is extrapolated to larger 
kR. The oscillating component of the cr± data is believed to 
be a numerical artifact. 



Lorentz force. For large kR, ka and — pf<Jl increase 
linearly with kR while fecii saturates to a constant value. 

The form of a and an can be understood in terms 
of a simple Born-limit calculation ^^|. The 2D Fourier 
transform of 1/r yields (approx) l/q and squares to yield 
da I dip ~ k/q 2 where q 2 — |k — k'| 2 = 2fc 2 (l — cosy). 
Since 1/(1 — cosy) diverges for small angles, angular inte- 
gration is dominated by the small angle cutoff at q w 1 jR 
and yields a ~ R. By constrast, since the extra factor of 
(1 — cosy) in the definition of cry precisely cancels this 
divergence, we obtain au ~ 1/fc. A nonzero oj_, however, 
can only be obtained by going beyond the Born limit. 

We can now use the fits to these numerical results as 
the input to Eqs. Q and J2J. The other input, A = 
K e (T)/aoT, is obtained empirically from the measured 
zero-field thermal conductivity in YBCO (extrapolated 
for the lowest T) 01- This quantity is a T-independent 
constant for T < 30 K but decreases for larger T where 
inelastic scattering becomes significant. Our calculated 
thermal conductivities are plotted (in scaling form) in 
Fig. [21 In each plot we show 15 curves for T ranging from 
15 K to 70 K. The low T curves (for which A ss const) 
satisfy scaling and therefore lie nearly on top of each 
other. At higher T, the curves deviate from scaling, pre- 
sumably due to the onset of inelastic scattering. Both 
the functional form of the scaling curves and the man- 
ner in which scaling is violated agree qualitatively with 
the mixed state thermal conductivity data measured in 
YBCO by Ong and co-workers Q. 

The form of our scaling curves can be understood as 
follows. In the small H regime, the mean free path is 
dominated by impurity scattering (^o <C £ v )- Since x = 
1/kR = \f~H j^T is small, the cross sections take on their 
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FIG. 2: Calculated longitudinal thermal conductivity (left) 
and thermal Hall conductivity (right) plotted in scaling form 
(in units defined in the text). The inset is the measured ther- 
mal Hall conductivity (in SI units) obtained by Ong and co- 
workers for YBa2Cu306.99 and reproduced here from Ref. 0. 



simple large-fci? form, f^(l/x) ss C|| and fj_(l/x) « c±/x, 
where C|| and cj_ are constants. Therefore, the scaling 
functions are F xx (x) — A and F xy (x) = cj_A x/tt and 
we find that the weak-field thermal Hall conductivity is 



(H, T) = C TVH 



(3) 



where Co = aokBC±A 2 /2irjEp. (For reasonable param- 
eter values, this yields C = 7.2 x Hr 3 W/mK 2 \/Tesla 
which is only a factor of 2 smaller than the experimental 
value reported in Ref.0.) As H grows, vortex scattering 
increases, increasing the skew scattering while reducing 
£ (and therefore n xx ). Near the point where £q ~ £ Vl 
the competition between the decreasing mean free path 
and the increasing skew scattering results in a peak for 
the K xy scaling curve. If the impurity scattering is suf- 
ficiently weak, this peak can occur for small enough H 
that the cross sections still take their simple large-fci? 
forms. Thus, for sufficiently clean samples, at temper- 
atures where scaling is satisfied, it follows from Eq. (0) 
that the peak height and location should be proportional 
to A 3 ! 2 and A^ 1 / 2 respectively. Recall, as well, that 
Co ~ A 2 . These are predictions that can be checked ex- 
perimentally by comparing results for samples of differing 
purity. Once vortex scattering dominates the mean free 
path (£ ~> i v )t K x V decreases with increasing H. As we 
push this model toward the strong-field (large x) regime, 
the transport cross section decreases from its large-fci? 
value, increasing £ and causing a leveling out and even- 
tual upturn in the scaling curves for both k xx and n xy . 
However, in this strong-field regime (Vll > jT), there is 
a magnetic field contribution to the quasiparticle density 
of states and so our picture of thermally-excited quasi- 
particles scattering from dilute vortices ceases to be valid. 



Note in particular, that the simple, yet previously un- 
expected, form of the weak-field thermal Hall conductiv- 
ity, n xy ~ T\/H, is now easily understood. Summarizing 
the preceding discussion, we write n xx ~ C v £ and K xy ~ 
K xx n v a±£ ~ C v £ 2 n v a±. In the weak-field limit, since 
quasiparticles are thermally excited, C v ~ T 2 . For small 
magnetic fields, the mean free path is dominated by im- 
purity scattering. Therefore £ ~ £o ~ n e (T)/C v ~ 1/T. 
The vortex density is just proportional to the magnetic 
field, n v ~ H. Our numerics yield a±_ ~ (k/pp)R ~ 
T I \/~H, which says that, aside from being small by a fac- 
tor of k/pp, the skew cross section is just proportional 
to the effective vortex radius, R. Putting it all together 
yields n xy ~ {T 2 )(l/T) 2 (H)(T/\/H) ~ T\/H. The sim- 
ple form of this result is due to the simple nature of the 
weak-field regime. Here we have thermally excited quasi- 
particles with a mean free path due to impurity scatter- 
ing. The only effect of the vortices is to add a small skew 
component to the scattering. The unusual \fH Hall re- 
sponse results because, while the number of vortices goes 
like H, the skew cross section per vortex goes like l/y/H. 

While our analysis appears to yield the correct func- 
tional form of the thermal conductivity, we do not expect 
our results to be quantitative due to the approximations 
we made in calculating the single vortex cross section. 
To obtain more quantitative results, one must include 
the Berry phase contribution, consider anisotropic Dirac 
nodes, and account (via empirical input) for the details 
of the vortex core. This is left for future work. 

We are very grateful to Y. Zhang and N. P. Ong for 
sending us their unpublished n xx data and allowing us to 
reproduce their n xy plot in Fig. [3 We thank D. Huse, A. 
Millis, T. Senthil, and J. Ye for helpful discussions. This 
work was supported by NSF Grant No. DMR-0201069. 
A. V. was supported by an MIT Pappalardo Fellowship. 



[1 
[2 
[3 
[4 

[5 

[6 
[7; 
[S 
[9 

[10 

[11 
[12 

[13 
[14 



K. Krishana et ai, Phys. Rev. Lett. 82, 5108 (1999) 

N. P. Ong et al, |cond-m at/9904160 (unpublished) 

R. Ocana et ai, J. Low Temp. Phys. 123, 181 (2001) 

Y. Zhang et ai, Phys. Rev. Lett. 86, 890 (2001) 

S. H. Simon and P. A. Lee, Phys. Rev. Lett. 78, 1548 

(1997) 

P. A. Lee, Phys. Rev. Lett. 71, 1887 (1993) 

M. J. Graf et ai, Phys. Rev. B 53, 15147 (1996) 

A. C. Durst and P. A. Lee, Phys. Rev. B 62, 1270 (2000) 

G. E. Volovik, Pis'ma Zh. Eksp. Teor. Fiz. 58, 457 (1993) 

[JETP Lett. 58, 469 (1993)] 

N. B. Kopnin and G. E. Volovik, Pis'ma Zh. Eksp. Teor. 

Fiz. 64, 641 (1996) [JETP Lett. 64, 690 (1996)] 

M. Franz, Phys. Rev. Lett. 82, 1760 (1999) 

I. Vekhter and A. Houghton, Phys. Rev. Lett. 83, 4626 

(1999) 

J. Ye, Phys. Rev. Lett. 86, 316 (2001) 

A. Vishwanath, Phys. Rev. Lett. 87, 217004 (2001); 

Phys. Rev. B 66, 64504 (2002) 



[15] O. Vafek et al., Phys. Rev. B 64, 224508 (2001) 
[16] Note that for an energy-dependent £, this simple relation 
between k xx and C v is only approximate. We verified 
that while the more exact relation (see Ref. llSfl modifies 
details, it preserves the qualitative form of our results. 



[17] Y. Zhang and N. P. Ong (private communication) 

[18] A. C. Durst et al. (in preparation) 

[19] R. Cleary, Phys. Rev. 175, 587 (1968) 

[20] M. Chiao et al, Phys. Rev. B 62, 3554 (2000) 



